ON AFFINE RIGIDITY 



STEVEN J. GORTLER, CRAIG GOTSMAN, LIGANG LIU, AND DYLAN P. THURSTON 

Abstract. We study the properties of affine rigidity of a hypergraph and prove a variety of 
fundamental results. First, we show that affine rigidity is a generic property (i.e., depends 
only on the hypergraph, not the particular embedding). Then we prove that a graph is 
generically neighborhood affinely rigid in d-dimensional space if it is ((i+l)-vertex-connected. 
We also show neighborhood affine rigidity of a graph implies universal rigidity of its squared 
graph. Our results, and affine rigidity more generally, have natural applications in point 
registration and localization, as well as connections to manifold learning. 



1. Introduction 

Suppose one has a number of overlapping "scans" of a set of points in some space, and 
that the corresponding points shared between scans have been identified. One naturally may 
want to register these scans and merge them together into a single configuration [23, inter 
alia]. Such a merging problem is called a realization problem. The study of the uniqueness 
of the solutions to such realization problems is known as rigidity. 

We model the combinatorics of this problem using a hypergraph B, with vertices repre- 
senting the points, and hyperedges representing the sets of points in each scan. The geometry 
of the problem is modeled with a configuration p, associating each vertex with a point in 
space. 

One natural setting is the Euclidean setting, where the scans are known to be related by 
Euclidean transforms. In this case it is sufficient to study just the case of a graph, where we 
think of each edge as its own scan with only 2 points. Unfortunately, many of the Euclidean 
problems are NP-HARD [18]. In this paper, we study what happens when one relaxes the 
problem to the affine setting, that is, one assumes that the scans are known to be related by 
affine transforms. Under this relaxation, much of the analysis reduces to linear algebra, and 
uniqueness questions reduce to rank calculations. We prove a variety of fundamental results 
about this type of rigidity and also place it in the context of other rigidity classes such as 
global rigidity and universal rigidity. 

We also specifically investigate the case of hypergraphs G that arises by starting with an 
input graph F, and considering each one-ring in F as a hyperedge in G. We call such a 
hypergraph the neighborhood hypergraph of F. Such neighborhood hypergraphs naturally 
arise when studying molecules [13] , when applying a divide and conquer approach to sensor 
network localization [20] and in machine learning [16]. 
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1.1. Summary of Results. We start by describing how affine rigidity in is fully char- 
acterized by the kernel size of one of its associated "affinity matrices". (This result was ffist 
shown by Zha and Zhang [26] . ) We show how this implies a number of interesting corollaries 
including the fact that affine rigidity is a generic property. That is, given a hypergraph G 
and dimension d, either all generic embeddings of G are affinely rigid in M'^ or all generic 
embeddings are affinely flexible in M'^. The speciflc geometric positions of the vertices are 
irrelevant to this property, as long as they are in sufficiently general position. Thus we call 
such a hypergraph generically affinely rigid in M'^. 

Next we relate affine rigidity in M"' to the related notion of universal Euclidean rigidity. 
A framework is universally Euclideanly rigid if it is rigid (in the Euclidean sense) in any 
dimension. In this context, we prove that affine rigidity in M.'^ implies universal Euclidean 
rigidity. 

We then prove the following sufficiency result: if a graph F is d + 1 (vertex) connected, 
then its neighborhood hypergraph is generically affinely rigid in W^; alternatively, we say 
that the graph F itself is generically neighborhood affinely rigid in M'^. In particular we will 
show that almost every non-symmetric equilibrium stress matrix for any generic embedding 
of F in M'^ will have co-rank d + 1 (i.e., rank v — d — 1). 

Putting these two results together, we show that if a graph is (i + 1 connected, then any 
generic embedding of its square graph in R'^ is universally rigid. This result is interesting, 
as very few families of graphs have been proven to be generically universally rigid. 

We give examples showing that many of the implications proved in this paper do not 
reverse. 

Finally we discuss some of the motivating applications. 

The main properties of frameworks of graphs and their implications proven in this paper 
are summarized below. 
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is generically globally rigid in M*^ 
is d -|- 1 connected ([10]) 

generically has non-symmetric equilibrium stress matrix of rank 

V — d — 1 (Proposition 5.8) 
is generically neighborhood affine rigid in M'^ (Proposition 5.9) 
is generically neighborhood universally rigid in M"^ (Corollary 4.2) 
is generically neighborhood globally rigid in M*^ (by deflnition) 



2. A RIGIDITY zoo 

In this paper we will consider several different rigidity theories. They all flt in to a unifying 
framework, which we now explain. 

Most generally, rigidity questions (of any type) ask if all of the geometric information 
about a set of points is determined by information from small subsets. In the usual Euclidean 
rigidity problem, we measure the distances between pairs of points. However, in other cases 
it is not enough to consider pairs of points for the small subsets; as such, we need to consider 
hypergraphs rather than just graphs. 
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Figure 1. Left: A hypergraph with 6 vertices and 4 hyperedges. Each hy- 
peredge is represented by a dotted elhpse enclosing a set of vertices. The 
hyperedges are: red {1,2,6}, green {2,3,5}, purple {5,6}, blue {4,5}. Right: 
The body graph of the hypergraph shown in the left. 

Definition 2.1. A hypergraph O is a set of v vertices V(0) and h hyperedges S{Q), where 
S{Q) is a set of subsets of V(0). We will typically write just V or £, dropping the hyper- 
graph G from the notation. 

There are natural ways to pass from a hypergraph back and forth to a graph. 

Definition 2.2. Given a hypergraph G, define its body graph B{Q) as follows. For each 
vertex in G, we have a vertex in -B(G). For each hyperedge h in Q and each pair of vertices 
in h we have an edge in -B(G). 

See Figure 1 for an example. 

Definition 2.3. Given a graph F, define its neighborhood hypergraph, written as A^(F) as 
follows. For each vertex in F, we have an associated vertex in A^(F). For each vertex in F 
we have a hyperedge in A^(F) consisting of that vertex and its neighbors in F. 

Definition 2.4. Given a graph F, its squared graph F^ is obtained by adding to F an edge 
between two vertices i and j if i and j share some neighbor vertex k. 

Lemma 2.5. For any graph T, B{N{r)) = F^. 

Proof. Immediate from the definitions. (See Figure 2 for an example.) □ 

Definition 2.6. A k-hypergraph G is a hypergraph where each hyperedge has exactly k 
vertices. For any A; € N and hypergraph G, let Bk{Q) be the fc-hypergraph whose hyperedges 
are all the subsets S of size k of vertices that are contained together in at least one hyperedge 
of G: 

S{Bk{e)) = {S\3he S{Q), Sch,\S\ = k}. 
For a vertex set S, the complete k-hypergraph on S, written Ki^^S), is the /c-hypergraph 
whose hyperedges are all (j^^^ subsets of S of size k. 

For instance, a 2-hypergraph is a graph, and -B2(G) is just the body graph -B(G). 
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Figure 2. Left: A graph with 6 vertices and 7 edges. Middle: Its neigh- 
borhood hypergraph with 6 vertices and 6 hyperedges: red {1,2,6}, pur- 
ple {1,2,3,6}, light blue {2,3,5}, dark blue {4,5}, orange {3,4,5,6}, green 
{1, 2, 5, 6}. Right: The body graph of the hypergraph in the middle. It is also 
the squared graph of the graph in the left. 



Definition 2.7. A configuration of the vertices V(0) of a hypergraph in M'^ is a mapping p 
from V(6) to W^. Let C^{V) be the space of configurations in R.'^. For p G C"^(V) and 
u G V(6), we write p{u) G M'^ for the image of u under p. 

A framework p = {p, 0) of a hypergraph is the pair of a hypergraph and a configuration of 
its vertices. C"^(G) is the space of frameworks {p, G) with hypergraph B and p G C"^(V(0)). 
We may also write p{u) for p{u) where p = {p, B) is a framework of the configuration p. 

A framework of a hypergraph has also been called a body-and-multipin framework [13]. 

Definition 2.8. Let G be a group acting on M'', for instance the group Eucl((i) of all 
Euclidean isometrics of M"^. Two frameworks p, cr G C'^(G) are G-equivalent if for each 
hyperedge h G ^^(6), the positions of the vertices in p can be mapped to their positions in a 
by an element of G depending only on h. That is, for each h & S, there is a, gh & G so that 
for each -u G /i, we have gh{p{u)) = a{u). The two are G-congruent if the positions of all the 
vertices in p can be mapped to their positions in a by a single element g & G (not depending 
on h). 

A framework p G C"^(0) is globally G-rigid if for any other framework a G C"^(G) that is 
G-equivalent to p, we also have that p and a are G-congruent. Otherwise we say that p is 
globally G-flexible in W^. 

Similarly, a framework p G G^(0) is locally G-rigid in if there is a small neighborhood U 
of p in G'^(9) so that for any a E U that is G-equivalent to p, we also have that p and a are 
G-congruent. Otherwise we say that p is locally G-flexible in M*^. 

A related notion of group based rigidity has been explored in the computer aided design 
literature [19]. 

In this paper, we are mainly concerned with the cases when G is either Eucl((i) or AS{d), 
the group of all invertible affine linear maps of M*^, in which case we speak about Euclidean 
or affine rigidity, respectively. But there are other interesting possibilities, like the group of 
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projective transformations. Another interesting case is when G is the group of dilations and 
translations of (with no rotations); this gives the theory of parallel-line redrawings [25]. 

In this terminology, Euclidean rigidity is the default: if the group G is not specified, it is 
the Euclidean group. In much of the rigidity literature, local rigidity is the default, and the 
qualifier "local" is dropped. However, in this paper this distinction is important and we will 
write "local" or "global" when the distinction is meaningful. 

Lemma 2.9. A framework (p, O) is locally (resp. globally) Euclideanly rigid iff the body 
framework {p,B{Q)) is locally (resp. globally) Euclideanly rigid. 

Proof. This easily follows from the fact that, for each hyperedge h G S{Q), the complete 
graph on vertices is globally rigid. □ 

Thus we only need to consider Euclidean rigidity for frameworks of graphs, not hyper- 
graphs. 

In the next section (Corollary 3.5) we will see that a framework is locally affinely rigid iff 
it is globally affinely rigid. Thus we can drop the local/global distinction for affine rigidity. 

Definition 2.10. For d < d', we view C^{V) as contained in C"^'(V) by the inclusion of M*^ 
as the first d dimensions of M'^ . 

Let G be a family of groups Gd acting on M'^, so that for d < d', Gd is the subgroup of Gd' 
that fixes M'' as a subset of M*^ . A framework p e C"^(0) is universally locally (resp. globally) 
G-rigid if it is locally (resp. globally) G-rigid as a framework in (0) for all d' > d. 

Note that universal rigidity of any sort implies rigidity of the same sort. 

Lemma 2.11. A framework p G C"^(0) is universally globally Euclideanly rigid iff it is 
universally locally Euclideanly rigid. 

Proof. For any two equivalent frameworks p in C"^(0) and p' in C"^'(0), Bezdek and Con- 
nelly [3] exhibited an explicit fiex between p and p' in C^+^'{Q). In particular, if p IS a 
rf-dimensional framework with a equivalent but non-congruent framework in d' dimensions, 
then this fiex shows that p is not locally rigid in ]R'^+°' . □ 

Thus we can also drop the local/global distinction in the case of universal Euclidean 
rigidity. 

Definition 2.12. A framework (p, F) of the graph F is neighborhood rigid (of any of the 
sorts above) if the corresponding framework (p, A^(F)) of the neighborhood hypergraph is 
rigid (of the same sort). 

For instance. Lemmas 2.9 and 2.5 tell us that neighborhood Euclidean rigidity of (p, F) is 
equivalent to the Euclidean rigidity of (p, F^). 

Definition 2.13. A configuration p in C^iV) is generic if the coordinates do not satisfy 
any non-zero algebraic equation with rational coefficients. A framework is generic if its 
configuration is generic. 

A property is generic in M*^ if, for every (hyper)graph, either all generic frameworks in 
C"^(0) have the property or none do. For instance, local and global Euclidean rigidity in M*^ 
are both generic properties of graphs and therefore for hypergraphs as well [2,8]. For any 
property P (generic or not) of frameworks, a (hyper)graph is generically P in if every 
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generic framework in C"^(G) has P. (For a non-generic property like universal Euclidean 
rigidity, there are (hyper)graphs that are neither generically P or generically not P.) 

The above adjectives are to be read in the reverse of the order we have defined them. 
Thus, for any framework, we may talk about 

(generic/0) (universal/0) (local/global) (Euclidean/affine) rigidity. 

where by we mean that this term has been dropped. 

3. Affine Rigidity in W'- 

We now move on the main focus of this work, affine rigidity, as defined in the previous 
section. Though the definitions start from a different point of view, this notion of affine 
rigidity is identical to the one defined by Zha and Zhang [26] and the concept is also informally 
mentioned by Brand [4]. Additionally, Theorem 1 below is essentially equivalent to [26, 
Theorem 5.2]. Our contribution here, described by the corollaries, is showing how affine 
rigidity fits in to the general scheme of rigidity problems. 

Lemma 3.1. Any general position framework of a complete {d + 2)-hypergraph is affinely 
rigid. 

Proof. Clear. □ 

Proposition 3.2. A framework (p, 0) in general position is affinely locally (resp. globally) 
rigid iff the associated framework (p, 5(^+2(0)) is affinely locally (resp. globally) rigid. 

(Compare Lemma 2.9.) 

Proof. First consider a hyperedge with less than d + 2 vertices. Any two generic embeddings 
in M*^ of these vertices are related by a (i-dimensional affine transform. Therefore this hy- 
peredge does not affect affine equivalence of these vertices, and so this hyperedge may be 
dropped without affecting affine rigidity. Next consider a hyperedge with k vertices, with 
k > d + 2. By Lemma 3.1, one can replace this hyperedge with (^^2) hyperedges correspond- 
ing to all subsets of c? + 2 vertices. The framework of the new hypergraph will be affinely 
rigid in M'^ iff the original one is. □ 

Definition 3.3. An affinity matrix of a framework {p, G) in C"^(G) is any matrix where each 
row encodes some affine relation between the coordinates of the vertices in a hyperedge of 
[p, G) as a homogeneous linear equation. That is, there are v columns in the matrix, the only 
non-zero entries in a row correspond to vertices in some hyperedge, the sum of the entries 
in a row is 0, and each of the coordinates of p, thought of as a vector of length v, is in the 
kernel of the matrix. 

An affinity matrix is strong if it encodes all of the affinely independent relations in every 
hyperedge of (p, G). 

Lemma 3.4. // the framework {q, G) is affinely equivalent to {p, G) then the coordinates of 
q are in the kernel of any affinity matrix for {p, G). Additionally, the converse is true if the 
affinity matrix is strong. 

Proof. Clear from the definitions. □ 

The kernel of an affinity matrix of a framework {p, G) G C"^(G) always contains the 
subspace of W spanned by the coordinates of p along each axis and the vector 1 of all I's. 
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This corresponds to the fact that any affine image of p is equivalent to p. If j9 is a proper d- 
dimensional configuration (with full d-dimensional affine span), these vectors are independent 
and span a (c? + l)-dimensional space. In particular, a generic framework of a hypergraph 
with at least d + 1 vertices in is proper, so for such frameworks the corank of any of its 
affinity matrices must be no less than d+1. 

The rank of strong affinity matrices fully characterize affine rigidity. 

Theorem 1. Let Q he a hypergraph with at least d + 1 vertices. Let (p, 0) he any proper, 
d- dimensional framework and let M he any strong affinity matrix for {p, 6). Then {p, 9) is 
affinely rigid in iff dim(ker(M)) = d + 1. 

Proof. By Lemma 3.4, for any other configuration q in C"^(V) such that (g, 0) is affinely 
equivalent to [p, B), the coordinates of q must be in the kernel of M. When dim(ker(M)) = 
d+1, the kernel of M contains only one-dimensional projections of p and the all-ones vector. 
Thus any (g, 9) that is affinely equivalent to {p, 0) must in fact be affinely congruent to 

(p,e). 

Conversely, if the corank is higher, then the kernel must contain an "extra" vector that 
is not a one-dimensional projection of p. Adding any amount of this vector to one of the 
coordinates of p must, by Lemma 3.4, produce a q such that (g, 0) is affinely equivalent to 
(p, 0) but not congruent to it. □ 

It is easy now to prove the following corollaries. 

Corollary 3.5. // {p, 0) is affinely glohally flexible in then it is affinely locally flexible. 

Proof. From the proof of Theorem 1, when {p, 0) is affinely globally flexible in M*^ there is 
an extra vector S in the kernel of a strong affinity matrix, and we can add any multiple of S 
to one of the coordinates of p to get an affinely equivalent but not congruent framework. □ 

Remark 3.6. In fact, for any two affinely equivalent frameworks (p, 0) and (g, 0), there 
is a continuous path of affinely equivalent frameworks in C"^(0) connecting them, namely 
((l-t)p + tg,0). 

Corollary 3.7. A framework (j9, 0) G C"^(0) is affinely rigid in M'^ iff it is affinely rigid 
when considered as a (degenerate) framework in M"^ for d' > d. 

Proof. Follows from the proof of Theorem 1. □ 

Thus there is no distinct notion of "universal" affine rigidity. 

Corollary 3.8. Affine rigidity in M'^ is a generic property of a hypergraph. 

Proof. The condition that M is an affinity matrix for (p, 0) is linear in the entries in M. 
The corollary then follows from Proposition A.l. □ 

Remark 3.9. Though we will not pursue the details here, one can use the concept of an 
affinity matrix to derive an efficient randomized discrete algorithm for testing generic affine 
rigidity of a hypergraph in M'^. To do this, one needs to use integers of sufficiently many 
bits, and do the arithmetic modulo a suitably large prime. The details parallel those in the 
global rigidity case [8, Section 5]. 
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4. Universal Euclidean Rigidity 

We now turn to universal Euclidean rigidity. To begin, we need the following technical 
definition: 

Definition 4.1. We say that the edge directions of a graph framework {p,T) G C^iT) are 
on a conic at infinity if there exists a symmetric d-hj-d matrix Q such that for all edges 
{u, v) of r, we have 

[p{u) -p{v)]'Q[p{u) ~p{v)] = 0. 

The edge directions of {p, T) are on a conic at infinity iff there is a continuous family of 
d-dimensional non- Euclidean affine transforms that preserve all of the edge lengths [6] . This 
is a very degenerate situation which is very easy to rule out. For example, if in a hypergraph 
framework [p, 6) some hyperedge in B has vertices whose positions in p affinely span M'^, 
then the edge directions {p,B{Q)) cannot be on a conic at infinity. 

Theorem 2. // a framework {p, 0) of a hypergraph B with p G C"^(V) is affinely rigid in M'^, 
and the edge directions of {p, B{Q)) are not on a conic at infinity, then (p, O) is universally 
Euclidean rigid. 

Proof. Let q G C"^'(V) be a configuration with d' > d such that (g, 6) is Euclidean equivalent 
in M'^ to 0). Then {q,Q) is affinely equivalent in M'' to (p, O). Since {p,Q) is affinely 
rigid in M.'^, from Corollary 3.7, we have that p is affinely congruent to q in M*^ and the affine 
span of q must be of the same dimension as that of p. Let R{q) be a rotation of q down 
to W^. Then R{q) is affine congruent in M.'^ to p and {R{q), O) is Euclidean equivalent in M*^ 

to (p,e). 

Since the edge directions of {p,B{Q)) are not on a conic at infinity, any ci-dimensional 
transform that is affine but not Euclidean must change some of the inter-vertex lengths 
within some hyperedge of 6. Putting this together with the fact that {R{q),9) and (p, 0) 
are affinely congruent and Euclidean equivalent, we can conclude that they are Euclidean 
congruent. Likewise q must be Euclidean congruent to p in M"^ . Thus we conclude that 
{p, O) is universally rigid. □ 

Corollary 4.2. Let Q be a hypergraph with at least d + 2 vertices. If a generic framework 
{p, 0) of a hypergraph O with p G C"^(V) is affinely rigid in M'' then {p, B) is universally 
rigid. 

Proof. Any generic framework of a hypergraph B with at least d + 2 vertices that is affinely 
rigid in M*^ must have at least one hyperedge h with at least d + 2 vertices, whose vertices 
in p have a ci-dimensional affine span. Thus {p, B{Q)) must include a generic framework of a 
{d + l)-simplex and thus cannot have edge directions at a conic at infinity. Then Theorem 2 
applies. □ 

There can be frameworks that are universally rigid but not affinely rigid in W^. (See 
Figure 3.) 

5. Neighborhood affine rigidity 

In this section we prove the following theorem about the generic neighborhood affine 
rigidity of a graph. 
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Figure 3. The framework in of the hypergraph on the left is not affinely 
rigid as each hyperedge (shown as a dashed elhpse) has only 3 vertices. But 
this framework is universally Euclidean rigid, as its body graph (right) is a 
Cauchy polygon. 

Theorem 3. Let T be a graph with at least d+1 vertices. IfT is {d + 1) -vertex- connected, 
then r is generically neighborhood affinely rigid in . 

We will build up the proof of Theorem 3 with a series of lemmas. 

Definition 5.1. An equilibrium stress matrix of a framework (p, F) of a graph in C"^(F) is a 
matrix Vt indexed by V x V so that 

(1) for all M, w G V, we have Vl{u, w) = Q{w, u); 

(2) for all u,w eV with u ^ w and {u, w} ^ £, we have Q{u, w) = 0; 

(3) for all M G V, we have J^w^v ^{u, w) = 0; and 

(4) for all M G V, we have J2w(^v ^{u, w)p{w) = 0. 

A non-symmetric equilibrium stress matrix of a framework (p, F) is a matrix that satisfies 
properties (2)-(4) above. 

Observe first that an equilibrium stress matrix (symmetric or not) Q of {p, F) is an affinity 
matrix of {p, A^(F)). From the properties of affinity matrices, the kernel of Q always contains 
the subspace spanned by the coordinates of p along each axis and the vector 1 of all I's. 

Definition 5.2. We say that a framework of a graph in C'^(F) has the convex containment 
property if 

(1) the one-ring of each vertex has an affine span of dimension d, and 

(2) except for d + 1 exceptional vertices, each of the other vertices in the framework is 
contained in the strict d-dimensional convex hull of its neighbors. 

Lemma 5.3. Let T be a graph with at least d -\- 1 vertices. Suppose F is {d -\- l)-connected. 
Then there exists a generic framework (g, F) in C"^(F) with the convex containment property. 

Proof. Pick any d-\- 1 vertices to be exceptional. Constrain the exceptional vertices to fixed 
generic positions in M*^ (at the vertices of a simplex). Associate generic positive weights Uij 
with each undirected edge ij. Find the "rubber band" configuration consistent with the 
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constrained vertices and these weights. Namely, find a framework (r, F) so that each non- 
exceptional vertex is the weighted linear average of its neighbors: 

Yl - r{j)) = 0, 

ieiv(i) 

where N{i) are the neighbors of vertex i. This involves solving d systems of linear equations, 
one for each component of r. Note that the resulting configuration r may not be generic. 

From [15], we know that if F is (rf + l)-connected and the constraints on the exceptional 
vertices and the edge weights u are generic, then no set oi d + 1 vertices in r will be contained 
in a ((i — l)-dimensional affine plane, giving us the first condition. 

By construction, any non-exceptional vertex in (r, F) must be contained in the convex hull 
of its neighbors. Again, from [15], the convex containment must be strict. 

Finally, we perturb each vertex in M'' to obtain a generic configuration in g G C"^(V). 
By the first convex containment condition, the convex hull of the neighbors of a vertex has 
non-empty interior, so a sufficiently small perturbation will maintain both conditions. □ 

Definition 5.4. Suppose that (g, F) has the convex containment property and is a non- 
symmetric equilibrium stress matrix for {q,T). We call a row of fl non- exceptional if its 
corresponding vertex is in the strict (i-dimensional convex hull of its neighbors. 

Lemma 5.5. Let T be a graph with at least d + 1 vertices. Suppose F is a {d+ l)-connected 
graph, and we have a framework {q, F) in C"'(F) with the convex containment property. 
Then there is a non-symmetric equilibrium stress matrix Q of {q,T), such that for every 
non- exceptional row i, we have the following property: If there is an edge connecting vertex 
i and vertex j , then is positive. 

Proof. All vertices have -|- 1 or more neighbors. For each vertex i, we can therefore find 
"barycentric coordinates": non-zero edge weights Uij on the adjoining edges so that 

jeN{i) 

If i is a non-exceptional vertex, due to the convex containment property we can choose the 
Uij to be positive. We then choose Qij = Uij for i ^ j and flu = — J2j ^ij- D 

Remark 5.6. This lemma is false if we require the stress matrix to be symmetric, because 
this prevents us from choosing Uij and Uji independently. 

Lemma 5.7. Let T be a graph with at least d + 1 vertices. Suppose F is {d + l)-connected, 
and we have a framework {q, F) in C"^(F) with the convex containment property. Then there 
is a non-symmetric equilibrium stress matrix Q of {q, F) with co-rank d + 1. 

Proof. From Lemma 5.5 we find for (g, F) a non-symmetric equilibrium stress matrix Q with 
the desired positive entries. We now show that Q has the stated rank. 

First remove the d+1 rows and columns associated with the exceptional vertices to create 
a smaller matrix Q'. Due to the sign pattern assumed in Q, as well as property (3) of any 
equilibrium stress matrix, Q' must be weakly diagonally dominant. 

Let us call a vertex EN if it has an exceptional neighbor and refer to its corresponding row 
in Q' as EN. Any EN row must be strictly diagonally dominant (since at least one non-zero 
off-diagonal entry of Q have been removed from this row). 

Since all entries corresponding to edges are non-zero, the reducible components of Q' 
correspond to vertex subsets that remain connected after the exceptional vertices have been 
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Figure 4. This framework in is not 3-connected but does have a non- 
symmetric stress matrix of high rank. 

removed. Each reducible component of Q' includes such an EN row, thus Q' must be full 
rank. (See, e.g., [24, Theorem 1.21].) 

Since fl' has co-rank 0, the co-rank of fl must be at most d + 1. It is no less since any 
equilibrium stress matrix must have a (d-l- l)-dimensional kernel spanned by the coordinates 
of q and the all-ones vector. □ 

Proposition 5.8. Let T be a graph with at least d+1 vertices. Suppose T is (d+l)- connected, 
and p is generic in C"^(V). Then there is a non- symmetric equilibrium stress matrix Q of 
{p, r) with co-rank d+l. 

Proof. From Lemma 5.3, there must exist a generic framework (g, F) in C"^(F) that has 
the convex containment property. From Lemma 5.7, (g, F) must have a non-symmetric 
equilibrium stress matrix of co-rank d+l. Thus from Proposition A.l, any generic framework 
{p, F) must have such a matrix as well. □ 

See Figure 4 for an example showing that the converse of Proposition 5.8 is not true. Here 
we show an neighborhood affinely rigid framework in with a non-symmetric equilibrium 
stress matrix of co-rank d+l = 3. This framework is not 3-connected. 

Note that from the proof of Proposition A.l it is clear that if F is ((i -|- l)-connected, then 
almost every non-symmetric stress matrix for almost every (p, F) in C"^(F) will have co-rank 
d + l. Moreover, each row of such a non-symmetric stress matrix of p can be constructed 
independently from the other rows, and we still expect to find this minimal co-rank. 

Proposition 5.9. Let T be a graph with at least d+l vertices. Suppose {p, F), a framework 
in C"^(F), has a non-symmetric equilibrium stress matrix fl that has co-rank d + l. Then 
(p, F) is neighborhood affinely rigid in M*^. 

Proof. Q is a, (not strong) affinity matrix of {p,N{r)) and so the proof follows that of the 
first direction of Theorem 1. □ 

Proof of Theorem 3. The theorem now follows directly from Propositions 5.8 and 5.9. □ 

See Figure 5 for an example showing that the converse of Proposition 5.9 is not true. Here 
we show a neighborhood affinely rigid framework in that does not have a non-symmetric 
equilibrium stress matrix of co-rank d+l = ?>. 

Remark 5.10. Generic global rigidity of a graph F in E*^ can be characterized either using 
the dimension of the kernel of a single symmetric stress matrix of a generic framework (p, F) 
or using the dimension of the shared symmetric stress kernel of a generic p: the intersection 
of the kernels of all stress matrices of p [8] . 
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Figure 5. This framework in does not have a non-symmetric equihbrium 
stress matrix of co-rank d + 1 = 3, but is (trivially) neighborhood affinely rigid. 

By contrast, the analogous statement is not true in the affine rigidity case. By "vertically 
concatenating" a sufficient number of non-symmetric equilibrium stress matrices of T), 
we can create a strong affinity matrix for (p, N{r)). The kernel of the vertical concatenation 
will be the shared non-symmetric stress kernel of (p, F), and the dimension of this kernel 
characterizes affine rigidity. Since the converse of Proposition 5.9 is not true, we see that 
neighborhood affine rigidity cannot in general be characterized by the rank of one single (say, 
generic) non-symmetric equilibrium stress matrix for (p, F). 

Note that there is a different sufficiency condition for affine rigidity given by Zha and 
Zhang [26]. Their condition is complementary to our condition (neither strictly stronger or 
weaker), and (like trilateration [7]) is greedy in nature. Their condition on generic frameworks 
of a hypergraph requires that one can walk between any two hyperedges such that for each 
sequence of two hyperedges along the walk, hi and hj share at least d+1 vertices. When 
translated to a neighborhood hypergraph A^(F), it states that one can walk between any two 
vertices such that for each sequence of two vertices along the walk, the neighborhoods 
of these vertices share at least d+1 vertices in F. 

Figure 6 shows a graph which clearly fails Zha and Zhang's condition, but is 3-connected, 
showing that their condition does not imply Theorem 3. It is not hard to construct examples 
in the opposite direction, as well. 

From Corollary 4.2, we also have the following corollary of Theorem 3. 

Corollary 5.11. Let T be a graph with at least d+1 vertices. If F is {d + l)-connected, 
then any generic framework of F^ in M.'^ is universally rigid. 

Remark 5.12. This corollary can also be proven without reference to affine rigidity and 
Corollary 4.2. In particular. Proposition 5.8 guarantees a non symmetric maximal-rank 
stress Q for {p, F), and then Q^Q is a symmetric, positive semi-definite, maximal rank stress 
for (p, F^). Universal rigidity then follows by a theorem of Connelly [5]. (See also [9].) 

We wish to highlight this corollary since the only other known (to us) class of graphs 
that are universally rigid for all generic configurations in M'' are the (i-trilateration graphs. 
(A (i-trilateration graph is one that can be obtained from a complete graph by successively 
adding vertices, each connected to at least d+1 old vertices.) 

See Figure 6 for an example of a framework whose square is universally rigid by this 
corollary but is not a trilateration graph. 

We also mention that a theorem of a related nature, showing a relationship between the 
connectivity of a graph and global rigidity in the squared graph, has been described by 
Anderson et al. [1]. 
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Figure 6. A drawing of a hexagonal lattice on the torus. (The vertices on 
the top edge should be identified with those on the bottom and similarly with 
the left and right, as indicated by the arrows.) This graph is 3-connected but 
its neighborhood hypergraph does not satisfy the sufficiency condition of Zha 
and Zhang [26], and its squared graph is not a 2-trilateration graph. 

6. Applications 

6.1. Registration. There are many applications where one has multiple views of some 
underlying configuration, but it is not known how these views all fit together. We assume 
that these views share some points in common, and this correspondence is known. (Of 
course in practice, establishing such a correspondence could in itself be a very challenging 
problem.) For example, in computer vision, one may have multiple uncalibrated laser scans 
of overlapping parts of some three-dimensional object. 

In our setting we model all of the points as vertices, and each of the views as a hyperedge. 
The geometry of the vertices in each hyperedge h is given up to some unknown transform Th 
from a relevant class. The goal in registration then is to realize the entire hypergraph up to 
the relevant congruence class. 

AfRne case: Suppose we wish to realize a framework {p, G) where we are given as input 
the geometry of each hyperedge h up to an affine transform A^. Theorem 1 tells us that 
if {p, 0) is affine rigid, then we can compute the realization just using linear algebra. In 
particular, we can use the data for each hyperedge to build its associated rows in a strong 
affinity matrix. Then we can solve for its kernel, giving us our answer p. 

If our hypergraph G happens to be the neighborhood graph of an underlying graph F, then 
one could also construct a (smaller) non-symmetric equilibrium stress matrix Q for (p, F). 
This is not guaranteed to work; even when {p, F) is neighborhood affinely rigid in M*^, the 
matrix Q may have co-rank larger than d + 1. But Theorem 3 states that if F is (c? + 1)- 
connected, this method will indeed work for almost every p in M'^ (and, in fact, using almost 
any non-symmetric equilibrium stress matrix for {p, F)). 

Euclidean case: The Euclidean framework registration problem is perhaps more natural 
and common. 
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When {p, G) is globally rigid in M'^, this problem is well posed, but it is general hard 
to solve, as the graph case includes the graph realization problem which is strongly NP- 
HARD [18]. 

When {p, B) is, in fact, also universally rigid there is an efficient algorithm: we can solve 
the Euclidean registration problem using semi-definite programming. One simply sets up 
the program that looks for the Gram matrix of an embedding of the vertices in W" (a semi- 
definite constraint on a Gram matrix) subject to the length constraints (linear constraints 
on the Gram matrix) [14]. Due to universal rigidity, one does not need to explicitly enforce 
the (non-convex) constraint that the embedding have a d-dimensional affine span [22]. 

When {p, 6) is, furthermore, affinely rigid then we can solve the Euclidean registration 
problem using linear algebra. We can simply reduce this problem to an affine registration 
problem above, and find p using the kernel vectors of an affinity matrix. This determines 
p up to some global affine transform. Moreover, for {p, 0) that is generically globally rigid 
in W^, we can solve a second (least squares) linear system to remove the unwanted global 
affine transform, leaving just the unknown global Euclidean transform (see Appendix B). 
This approach is morally the same affine relaxation used in the initialization step of the 
registration method of Krishnan et al. [12] (though in their case, they think of the inter- 
patch transforms as the unknown variables instead of the point positions). 

6.2. Global embeddings from edge lengths. Similar approaches have been applied to 
the (NP-HARD) problem of solving for the framework of a graph given its edge lengths. In 
these approaches one first attempts to find local d-dimensional embeddings for each one-ring 
of the framework up to an unknown local Euclidean transform. This step alone is NP-HARD 
and can fail. But assuming this step is (approximately) successful one can reduce the rest 
of the problem to the Euclidean registration problem above. 

In the As-Affine-As-Possible (AAAP) method [11,27], this was done using what is essen- 
tially a strong affinity matrix. In the Locally-Rigid-Embedding (LRE) method [20] this was 
done using a non-symmetric equilibrium stress matrix. Both approaches then removed the 
global affine transform using the least squares linear system described in Appendix B. 

6.3. Manifold learning. Many of the ideas of affine rigidity first appeared in the context of 
manifold learning. Suppose one has c?- dimensional smooth manifold which is a topological 
(i-ball embedded in a larger D-dimensional space M^. Also suppose that one has a set V of 
sample vertices on the manifold. In manifold learning, one first connects nearby samples to 
form a proximity graph F. One then looks for a framework [p, F) of this graph in M*^ that in 
some way that preserves some of the geometric relations of the points in M.^. This is used 
to represent a parametrization of Ai. 

To compute the coordinates p, the Locally Linear Embedding (LLE) method [16] builds a 
matrix Q with structure similar to a non-symmetric equilibrium stress matrix. In particular, 
row i encodes one affine relation between vertex i and its neighbors in M^. Then (after ignor- 
ing the all ones vector) the smallest d eigenvectors of Q^Q are used to form the coordinates 
of p in M*^. Unfortunately, since the original embedding is in M.^, for a graph of high enough 
valence and assuming no noise, Q must have a kernel of size at least D + 1, which is much 
larger than d + 1. Thus it is not clear how the numerically smallest d + 1 eigenvectors will 
behave. The paper suggests to add an additional regularization term, possibly to address 
this issue. 
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A follow up to the LLE paper [17] describes a PCA-LLE variant where a ci-dimensional 
local PCA is computed to "flatten" each one-ring before calculating its corresponding row 
in the matrix Q. Thus Q is designed to represent c?- dimensional affine relations between the 
points. 

The Local- Tangent-Space- Alignment (LTSA) method [28] is an interesting variant of PCA- 
LLE. In this method, a, v x v matrix is formed that is the Hessian of a quadratic energy. 
This energy vanishes only for vectors that are affinely equivalent to each of the flattened 
one-rings. Thus this matrix plays the role of a strong affinity matrix. It is in this context 
that Zha and Zhang investigated the rank of this matrix and affine rigidity [26]. 

In all of these methods, an understanding of affine rigidity is important. In particular it 
tells us what the rank of the computed matrix would be if the original d-dimensional manifold 
was in fact embedded in M*^. For example, if such a framework was not affinely rigid in M'^, 
then the kernel would be too big, and we would not expect a manifold learning technique 
to succeed. However, in manifold learning the embedding has an affine span greater than 
d and the analysis becomes more difficult. The kernel of a strong D-dimensional affinity 
matrix is too large, while the kernel of a strong affinity matrix for the locally flattened 
conflgurations contains only the all-ones vector but it is hoped that the numerically next 
smallest d eigenvectors are somehow geometrically meaningful. For an analysis along these 
lines, see [21]. 

Appendix A. Matrix rank 

For completeness, we recall the necessary material for determining the generic matrix rank. 
This material is standard. For a more detailed treatment, see, e.g., [8, Section 5]. 

We will consider the general setting where there is a set of linear constraints that must 
be satisfled by a vector m G M". The entries of m are then arranged in some flxed manner 
as the entries of a matrix M, whose rank we wish to understand. The linear constraints 
are described by a constraint matrix with n columns: C{p). Each of the entries of C{p) is 
deflned by some polynomial function, with coefficients in Q, of the coordinates of an input 
configuration p. We wish to study the behavior of the rank of M as one changes p. 

We apply this in the proof of Corollary 3.8 where the constraints C specify that the 
matrix M is an affinity matrix for {p, 0) and in the proof of Proposition 5.8 where the 
constraints C specify that the matrix M is a non-symmetric equilibrium stress matrix for 

(p,r). 

Proposition A.l. Suppose that for some generic p, there is a matrix M of rank s consistent 
with C{p)m = 0. 

Then for all generic p, there is some matrix M of rank > s consistent with C{p)m = 0. 
Proof. To prove this proposition we first need the following lemma. 

Lemma A. 2. Let M(7r) be a matrix whose entries are polynomial functions with ratio- 
nal coefficients in the variables tt G M". Let r be a rank achieved by some M(7ro). Then 
rank(M (tt)) > r for all points vr that are generic in M". 

Proof. The rank of the M(7r) is less than r iff the determinants of all of the r x r submatrices 
vanish. Let ttq G M" be a choice of parameters so M{no) has rank r. Then there is an 
r X r submatrix T('7ro) of M{7Tq) with non-zero determinant. Thus det(T(7r)) is a non-zero 
polynomial of tt. For any vr with rank(M(7r)) < r, this determinant must vanish. Thus, any 
such TT cannot be generic. □ 
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Next we recall that for a non-singular n x n matrix C, 

(1) a.di{C) = det{C)C-\ 

where adj C is the adjugate matrix of C, the conjugate of the cofactor matrix of C. In 
particular, adj C* is a polynomial in C. 

For a given p, let t{p) be the rank of C{p). Let t := maxpt(p). By Lemma A. 2 this 
maximum is obtained for generic p. 

For each p we add a set if of — t additional rows to C{p) to obtain a matrix C{p, H), 
and determine m by solving the linear system C{p,H)m = b where b G MJ^ is a vector of 
all zeros except for a single 1 in one of the positions of a row in H (if there are any rows 
in H). This m is then converted to a matrix M{p, H). M{p, H) is well-defined iff this linear 
system has a unique solution, i.e., iff C(p, H) has rank n. Note that this happens for generic 
p and H. 

Let Pq be generic and have a compatible matrix Mq with rank s, as in the hypotheses 
of the proposition. Find a set Hq of additional rows so that C{po,Hq) has rank n and 
C(po, HQ)mQ = b. Let C{p,H) be an n x n submatrix of C{p,H) so that C{po,Hq) is 
invertible. (C* necessarily uses t rows from C(p) and all rows of H.) Define b similarly, let 
7h{p,H) := adj(C*)6, and let M{p,H) be the associated matrix. 

By Lemma A. 2, the rank of C{p, H) is equal to its maximum value, n, at all points (p, H) 
that are not zeros of a polynomial Pi{p,H) := det C{p, H). Moreover, when Pi{p,H) ^ 0, 
the linear equation defining M{p, H) has a unique solution and the adjugate matrix M(p, H) 
is a scalar multiple of M(p, H). In particular we have assumed (po, -f^o) is not a zero of Pi 
and thus M{po,Ho) has rank s. By Lemma A. 2 again, the rank of M{p,H) is less than s 
only at the zeros of a non-zero polynomial P2{p, H). 

For any generic p, there must be some generic point (p, H). Such a generic (p, if) cannot 
be a zero of Pi or P2 and thus M(p, H) and M(p, if) must have rank no less than s. □ 

Appendix B. Removing the Affine Transform 

Suppose one has solved for g - a configuration in M*^ - up to an unknown global afiine 
transform A of the true configuration p: p = A{q). Given a set of edge lengths for p, it is 
possible to compute A up to an unknown global Euclidean transform. This approach was 
described by Singer [20]. 

In particular, let L be a c/x matrix representing the linear portion of A and let G := L^L. 
Now consider the following set of linear equations (in the unknowns of G): For each 

pair of vertices i,j, whose edge lengths are known, we require 

(2) {qi^) - q{3)YG{q{i) - g(j)) = {p{i) - p{3)Y{p{^ - p{j)) 

(Since we have more constraints than unknowns, for numerical purposes we would solve this 
as a least squares linear system in the unknown G.) 

The only remaining concern is whether this system has more than one solution. The 
solution to Eq. (2) will be unique as long as our set of edges with known lengths do not lie 
on a conic at infinity. 

Fortunately, we can conclude from Proposition 4.3 of [6] that if our known lengths form 
a graph B{Q) with minimal valence at least d and p is generic, then the edges do not lie on 
a conic at infinity. This property holds for any hypergraph that is generically globally rigid 
in E"^. 
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Using Cholesky decomposition on G then yields L up to a global Euclidean transform. 
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